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Applying the London theory we study curved vortices produced by an external current near and 
parallel to the surface of a type II superconductor. By minimizing the energy functional we find 
the contour describing the hard core of the flux line, and predict the threshold current for entrance 
of the first vortex. We assume that the vortex entrance is allowed due to surface defects, despite 
the Bean-Livingston barrier. Compared to the usual situation with a homogeneous magnetic field, 
the main effect of the present geometry is that larger magnetic fields can be applied locally before 
vortices enter the superconducting sample. It is argued that this effect can be further enhanced in 
anisotropic superconductors. 
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I. INTRODUCTION 

Surface barrier effects in type II superconductors have 
been predicted by Bean and Livingston 1 and de Gennes^ 
The entry of flux lines into a planar type II supercon- 
ductor situated in an external magnetic field H ext par- 
allel to its surface is opposed by a strong surface bar- 
rier when H ext = H c \, the first critical field. There- 
fore the entry of flux lines could occur at a field value 
Hext = Hs ~ H C 2 H c i, where H C 2 is the second criti- 
cal field. These surface barrier effects have been observed 
experimentally in the 60's on lead thallium alloys^ and on 
niobium metal,— and make it difficult to measure directly 
thermodynamic properties of the superconductor. 

Typically surface barriers are reduced due to surface 
disorder, which creates large local magnetic fields and 
allows for nucleation of vortices. Suppression of surface 
barriers for flux penetration was observed on YBaCuO 7 
and in BiSrCaCuO whiskers- due to heavy ion irradi- 
ation. In ellipsoid-shaped YBaCuO it has been argued 
that due to roughness of submicrometer order the surface 
barrier does not push the penetration field Hs above H c \ 
but only lowers the rate of vortex entry^ 

Another source for the delay of the entry of flux lines 
into superconductors is the "geometrical barrier" , 10 i n 
which is particularly important in thin films of constant 
thickness (i.e., rectangular cross section). This effect is 
absent only when the superconductor is of exactly ellip- 
soidal shape or is tapered like a wedge with a sharp edge 
where flux penetration is facilitated. The resulting ab- 
sence of hysteresis in wedge-shaped samples was nicely 
shown by Morozov et ak- 
in this paper we study another source for delay of en- 
trance of flux lines, due to inhomogeneity of the external 
magnetic field. In particular we consider magnetic field 
produced by an external current / flowing parallel to the 
surface of a type II superconductor, see Fig. 1. The 
magnetic field produced by the external current enters 
the sample as curved vortices at sufficiently large current. 
We find that the entrance of the first line occurs when the 
induced magnetic field at the surface at the position clos- 
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FIG. 1: Curved flux line near the surface of a superconductor 
enabled by an external current I. 



est to the wire already exceeds the bulk critical field H c \ . 
This delay in entrance of the curved vortices occurs due 
to geometrical reasons: The entry and outlet points are 

i 2 

associated with an energy cost ~ , where A is the pen- 
etration depth, (f>Q is the flux quantum, and /in is the free 
permeability. Note that ^q/mo^s = 0.2464 10 6 K — fxm, 
implying that in typical superconductors this is a large 
energy scale. In addition the spatially averaged mag- 
netic field experienced by the vortex is lower than the 
the maximal one occurring closest to the wire. Consid- 
ering those effects in an actual calculation we find how 
large a magnetic field can be applied locally without in- 
troducing vortices into the sample. 

This implies that application of magnetic field by an 
external current near the SC can be convenient for exper- 
iments demanding sizable magnetic fields in the vortex- 
free state. As such an experiment we mention the 
London-Hall effect^ Whereas this effect was observed 

it is now interesting to 
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in regular superconductors ,° 
measure it in high temperature superconductors. Typi- 
cally H c \ is quite low in these materials and therefore 
vortices penetrate the sample at very low homogeneous 
magnetic fields; hence our geometry can be useful. How- 
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ever other surface effects seem to be an additional ob- 
stacle for the observation of the London effect in high 
temperature superconductors.— 

A parameter which we leave out of consideration in 
this work is anisotropy of the superconductor, which is 
particularly important in high temperature layered su- 
perconductors. In the case of strong anisotropy addi- 
tional complications enter the problem even in the case 
of uniform magnetic field, where the direction of the vor- 
tices deviates from the direction of the external magnetic 
fieldi^i For certain (elliptical) treatment of the short dis- 
tance cutoff the vortices can have two different direc- 
tions, corresponding to two degenerate minima in the 
free energy^ 

We argue that strong anisotropy is expected to have 
important effects in our geometry, increasing further the 
maximal local magnetic field allowed before curved flux 
lines penetrate the sample. Consider the case where the 
c axis of a uniaxially anisotropic superconductor corre- 
sponds to the direction x in our geometry, c |j x. In this 
case the surface of the superconductor, parallel to the 
external wire, corresponds to an ab-plane. In the limit of 
strong anisotropy X a b *C A c the bulk critical field parallel 
to the surface H c \ = -. — ^ — r- log % becomes very small. 
On the other hand, the entry and outlet points of the flux 
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a single flux line. 



line are associated with a large energy cost 
pendent of A c . Therefore we expect the maximal surface 
magnetic field before the entry of the first vortex to in- 
crease relative to H c \ as a function of X c /X a b- We leave 
a detailed treatment of anisotropy in this geometry for a 
future work. 

The paper is organized as follows. In Sec.|TT]we formu- 
late the problem and obtain expressions for the magnetic 
field and free energy within London theory. In Sec. IIIII 
we present and discuss the numerical results for the min- 
imization of the free energy as a function of vortex con- 
tour. In Sec. HI! and Sec. lIII Al we consider the simpler but 
unrealistic case of a wire with zero width (i.e. < A). In 
Sec. lIII Bl we generalize to wires with finite width. Sec. lIVI 
contains conclusions. Some details about the derivation 
of the expression of the free energy are relegated to the 
appendix. 



II. FORMULATION 

Suppose that a type-II superconductor (SC) occupies 
the region x > and magnetic field is induced by an ex- 
ternal current / flowing along a wire of zero cross-section 
at {x, z) = (— d, 0), see Fig. ([lj. Our main object under 
consideration is a curved flux line lying in the plane y = 0. 
Let 7 denote the closed contour in Fig. @ consisting of 
the axial line of the flux line T and a line Ti symmetric 
to r with respect to the plane x = 0, corresponding to an 
image vortex. Upon further increasing the current a lat- 
tice of curved vortices is expected to form along the wire. 
However here we shall concentrate on small currents and 




FIG. 2: r is the axial line of the vortex line, 
contour 7 is V + Ti . 
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In the type II limit, where the coherence length £ is 
much shorter than the penetration depth A, the total 
free energy at zero temperature is given by 2 



F[T] = I 



/ d 3 r[H 2 + 6»(x)A 2 (V x H) 2 } 

Jr>f 



d 3 rA-j cxt 



(1) 



Here j C xt = —IS(x + d)8(z)y, I is the applied current 
through the wire, A is the vector potential H — W x A 
and 9(x) is the unit step function. The integral J r> ^ is 
carried out in all space outside of the vortex "hard core" 
T. We assume that the radius of curvature of T is larger 
than £ at any point in T. Note that at x = there is an 
apparent kink in 7, however this should be though of as 
a kink only for length scales large compared to £. 

The corresponding equations for the magnetic field H 
are the Maxwell equation, V x H = j C xt for x < 0, and the 
London equation, (1 - A 2 V 2 )i?(r) = ^ j v dr'5 3 {r- f) 

for x > 0. For all x we also have V-if = 0. In addition we 
impose appropriate boundary conditions at x = 0: the 
magnetic field is continuous, and no supercurrent flows 
perpendicular to the surface: j x = (V x H) x — 0. To 
construct a solution we use the functions 
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Here k 2 — k y y + k z z, k 2 — y^y + an d T (k) — 



Vfc 2 + A- 2 . For any A(fc 2 ), B(fe 2 ), the function H hom 
satisfies the homogeneous equations V x H hom = for 
x < 0, and (1 - \ 2 V 2 )H hom = for x > 0. The func- 
tion 7J 7 satisfies the London equation (1 — A 2 V 2 )_ff 7 (r) = 
^ J dr*5 3 (r — f*) in all space. The function Hp^* sat- 
isfies Maxwell equation V x Hii^>(r) = Joxt f° r Jext = 
/'<5(a; + d')8(z)y for all space. 

Defining the surface 2-dimensional Fourier transform 

H™ rf (k 2 ) = J dydze-^H^O, y, z), for fx = 7 , {/', a"}, 
one finds 



r ( k 2)\l~'J 



dr e 



7 2^ A 2 J 7 T(fc 2 ; 

= S(k y )wle-^(i sgn fc z ,0,sgn d). (3) 

The solution of the equations satisfying the desired 
boundary conditions is obtained by adding together the 
functions in Eq. ©, and solving for A(k 2 ) and B(k 2 ) 
to give continuity, ft is convenient to include an image 
current at x = — d. The total magnetic field is 

H = Ho + H v + H s , 

H = 9{-x){H I4 + H- It - d ) + H s o, H s0 = H^ m Bo , 
H v = 9(x)H„ H s = H\° m Bi , (4) 
where 

fsurf , 



A Q (k 2 ) = 



Aifa) = B 1 {k 2 ) = 



2[Hj^(k 2 )] z 
ik z (r(k 2 ) + k 2 

[H™ r f(k 2 )} x 



B (k 2 ) = -A (k 2 ) 



k 2 [k 2 + T(k 2 )}' 



r{k 2 y 

(5) 



In the absence of vortices the magnetic field is given by 
H a . It is plotted in Fig. © for d = 5 A. 




FIG. 3: Field lines of the vortex-free solution Ho(x,z) for 
d = 5A [direction of field lines correspond to anticlockwise 
rotation around (x,z) — (— d, 0)]. 



The total free energy as function of T is obtained by 
substituting the magnetic field Eq. (0| into the free en- 
ergy Eq. (P). We obtain 



^ J d\[H 2 a + 6(x)\ 2 (V x H Q ) 2 - 2A ■ j cxt ] , 

Jd 3 r[H 2 + 9(x)X 2 ^xH l ) 2 ] 7 i = v,s, 
-Ho / d 3 r(A v +A S ) • J cxt . (7) 



Here 
F 

Fi 

Fext 



Here Hi — V x Ai, i — 0, v, s. All mixed terms between 
Hq, H V1 H s vanish. For the vanishing of mixed terms in- 
volving H v see Ref. [|1 p. 5 79. We prove the vanishing of 
the remaining crossed terms between Ho and H s in the 
appendix. 

The term Fq is the energy of the system without vor- 
tices. To evaluate it we introduce a finite wire radius 
a«A,d, and assume the external current flows in a thin 
shell of this radius. Note that Fq scales linearly with the 
length of the wire, L y . The result of a calculation, using 
the methods of the appendix, is 



£o = Vol 2 

L. 



g(v) 



2tt 
dx 



log(2d/a)+g(d/X) 



-2x 



+ \J x 2 + y 2 



(8) 



We can infer from it the repulsive force per unit length 
Qfp 1 - < between the wire and the SC. It is plotted 

in Fig. © (for a/A = 0.01). Using g(y -> oo) = ±, 

g(y — > 0) = |logj^, we may identify two regimes, (i) 

d > A: Here g — ► and %^ _> _ J£o|L. i n agrement 



with Ampere force law, this corresponds to a repulsive 
force per unit length between two wires 2d apart carry- 
ing current / with opposite direction. This is the ori- 
gin of the levitation effect. The second wire corresponds 
to the term H-1-4 in the solution for H , see Eq. ((U); 
(ii) d -C A: The 1/d divergence in the force is cutoff by 
A. The limiting repulsion force per unit length as the 

wire approaches the surface is ^jf 2 - — > — ^x c o where 
c - 0.665. 

The term F ext accounts for the interaction between the 
vortex and the external current. Using j ext = —18{x + 
d)5{z)y we have 

/oo 
dy[A v (-d,y,0)+A s {-d,y,0)} v . (9) 
-00 

The contour of integration (2, y, z) — (0, —00, 0) — > 
(0, 00, 0) corresponds to the external current. Physically 
the wire should be closed into a loop, and we may close 
the contour of integration e.g. in the xy plane from the 
x — > —00 side. Then, using Green's theorem we obtain 



F = F Q + F V +F S + F e , 



(6) 



F ex t — Vol 



dx 



dy[H v (x,y,0) + H s (x,y,0)} 2 
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FIG. 4: Repulsive force between the superconductor and 

the wire in the absence of vortices. The dimensionless force FIG. 5: Contours corresponding to local minimum of F for 

-(2n\/L v ii I 2 )d d F behaves as A/(2d) for d > A, and goes d = ioa, and for the specified currents. We assume that V has 

to a constant c ~ 0.665 for d < A. the reflection symmetry z -> -z, and plot F only for z > 0. 



Note that H v (x,y,0) vanishes at x < 0. Using the for- 
mula for H Sl Eq. we obtain 



F, , , -— /" dr z / dke- kd cos(kr z ) x 



r Jo 



(10) 



We used the identity § dr ■ VJ-(r) = which holds for 
any continuous function J 7 if 7 is a closed contour. In 



r (fc)k,| 

The terms F v and F s have been derived in Rcfs. 



this calculation T(r) = e 1 r *sgn(r x ) (1 — e 



kg, 



2^0 - ^ 



dr,; / 



, exp(-|r-/|/A) 
,; 87rA 2 lf-f*| ' 



F s = ^ f dr z f ^yW(r-f). 
^0 Jr -M 



(11) 



The term F v is sensitive to the short distance cutoff £. To 
account for the cutoff we restrict the contour integration 
to \r — r*| > £. Here the anisotropic kernel for i 7 ^ is 



1 



2vrA 2 



dk 11- 



r(k) 



Jo(k\r z \), 



where Jq(x) is a Bessel function, and this integral can be 
done and expressed in terms of other Bessel functions. 14 
Note that U (s) (r* -> 0,r z -> 0) = (27rA 3 )-\ hence there 
is no need to regulate F s with a cutoff. 

Different than the usual case with a uniform magnetic 
field, in our problem the energy F — Fq + F v + F s + F ext 
is a function of the contour Y and is minimized for a par- 
ticular contour which we need to find. To this end we 
minimize F[T] numerically, approximating T by a poly- 
line having 2M equal length sides (M = 8 in most simu- 
lations). We assume that V has the reflection symmetry 
z — > — z. This leads to a M + 1-dimensional parameter 
space in which we search for the minimum of F. For an 
example see Fig. ([3]). In all our calculations £ = .001 A. 



III. SURFACE BARRIER 

We find that the free energy contains a surface energy 
barrier. From this section we shall disregard the vortex 
independent term of the free energy, F — > F v + F s +F ext . 
For later comparison we briefly discuss the case with ho- 
mogeneous magnetic fieldi-i Consider a semi-infinite type- 
II superconductor with a flux thread within it, parallel to 
the surface and to the external magnetic field H ext (|| z). 
The line energy / = F/L z (L z is the length of the vortex 
taken to be parallel to z) as function of the distance from 
the surface xq, is given bj*^ 

f(x ) = 0o [H ext e- X °' x - h(2x ) + H cl - H ext ] . (12) 

Here h(r) = 27r ^ A !i Ko(j) is the function giving the field 
at distance r of a single straight flux line, H c \ = \h{£) = 
; log h, and Kn is the zero-order Bessel function. 

471740 G £ 

The term (j)oH ext e~ x °/ x describes the interaction of the 
line with the external field and associated screening cur- 
rents. It is a repulsive term. The term —4>oh(2xo)/2 
represents the attraction between the line and its image. 
When Hq ~ H c \ there is a strong barrier opposing the 
entry of a line. We can understand this barrier as fol- 
lows: When H ext = H cl , f(x = 0) = f{x = 00) = 0. 
If we start from xq large and bring the line closer to the 
surface, the repulsive term (~ e~ x °/ x ) dominates the im- 
age term (~ e~ 2x °/ x ). Thus / becomes positive and we 
have a barrier. The barrier disappears, however, in high 
fields. When H > H$ = 0o/47rA£, the slope df /dx \ Xo ^^ 
becomes negative. H$ is of the order of the thermody- 
namic critical field H c i. The conclusion is that, at field 
H < Hs, the lines cannot enter in an ideal specimen (al- 
though their entry is thermodynamically allowed as soon 
as H > H c i). However this picture is modified in ex- 
periment due to surface inhomogeneities producing local 
large magnetic fields, and allowing vortices to enter the 
sample above H c \. 
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A. Results for wire with zero width 

Wc find a similar energy barrier for the entrance or exit 
of a curved vortex in our geometry with an external cur- 
rent rather than an homogeneous external magnetic field. 
This barrier can be visualized in the curves in Fig. ([6]) 
(except for the diamonds). Note that typically the bar- 
rier hight A is of order A ~ </>o/Mo ^ ^ T c , where T c is the 
critical temperature of the SC. This implies rather small 
tunneling probabilities e~ A / T -C 1 which prevents entry 
of vortices for clean surfaces. However for strong disor- 
der, vortices can enter more efficiently via nucleation at 
impurity sites. The contours corresponding to the mini- 
mum of the curves with stars and squares are plotted in 
Fig. ([5]). In all our calculations £ = .001A. 
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FIG. 6: Evolution of surface barrier as function of external 
current for d/X = 10. When I < I c0 (diamonds, I = 
the force on the line points always towards the surface. When 
Jco < I < Id (stars, / = 22^^) there exists a meta-stable 
minima with positive energy. When 7 c i < I < Is (squares 
and triangles I — 26, 80 x the minimum energy is neg- 

ative, but a barrier opposes the entry of the flux line. In 
each point in this plot we have minimized numerically F with 
respect to the contour V at fixed xo- 



Figure (|6|) implies the following picture. For infinites- 
imal current there is no stable vortex configuration. As 
the current increases we identify three threshold currents 
I c q < I c i < Is- When the current exceeds I c q a meta- 
stable minima with F > occurs. When the current 
exceeds I c ±, the minimum energy changes sign, F < 0, 
but still there is an energy barrier for the entrance of a 
flux line. When the current exceeds Is the barrier disap- 
pears. 

In Fig. |(7J) we investigated the dependence of I c q and 
I c i on the distance to the wire d. In the limit d X 
the results for I c \ are consistent with the formula I c \ — > 
ndH c i [see diagonal dashed line in Fig. ([7])]. The behavior 
of I c q in that limit shows that the region of metastabil- 
ity J c o < I < I c i is very narrow. This behavior appears 
in sharp contrast to the case of uniform magnetic field 
even in the limit d ^> A: We recall that Eq. fT2")) predicts 
metastable solutions for infinitesimal homogeneous mag- 
netic field H ext . These states live far from the surface as 



H ext becomes smaller. This effect is absent in our geom- 
etry both due to the fact that the effective external mag- 
netic field created by the wire decays at long distances 
from the surface and due to the line energy for pene- 
tration a long distance into the SC. In the other extreme 
limit d <C A we observed from the numerical solution that 
the contour 7 can be approximated by a circle centered 
at the origin. Making this assumption we can calculate 
jrircfa = 10 - 6981 ^ (x ~ 0.72A, Ffi X/ 4% = 0.1715), 

and I^ rcle = lO 749 ^ (x Q = 1.27A, F = 0) in the limit 
d — ► 0. This approximation is in reasonable agreement 
with the actual solution as the horizontal dashed lines 
show. 




FIG. 7: Dependence of threshold currents 7 c o and I c i on d/X. 



The shape of the contour changes as function of d. 
In Fig. ([8]) we plot the extension of the contour in the 
x and z directions. We fitted the numerical results for 
Xq with an empirical formula xq/\ = c + log(d/A) with 
c ~ 1, implying that the penetration of the vortex is of 
order A for all d. On the other hand it appears that zq 
grows linearly as function of d. In the limit d — > we 
have Xo/X — > 1.26 and zo/X — ► 1.43 (implying that the 
circular contour is only an approximation). 
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FIG. 8: Extensions of the curved flux line along x and z as 
function of d/X, at I = I c \. 



For disordered surfaces, the present geometry can be 
useful for application of large magnetic fields on a SC 
sample in a vortex free state. The maximal magnetic field 
that can be applied in a vortex free state using the wire 
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geometry is Ho[x = z = 0] [see Eq. ((4])] at current I&, 
In the limit d > A this magnetic field coincides with the 
bulk first critical field H c \ w 47r ^° A 2 log(A/g), however 
at smaller d the magnetic field at the surface increases. 
This is shown in Fig. @ where we plot the magnetic 
field H surface = [i? (0,0,0)] 2 = [i? s0 (0+,0,0)] z given in 
Eq. (J4j) at the current I c \, which we calculated above as 
function of d. Note that the field enhancement is small 
for d > 3 A (H surface = 2H cl for d = 3A). 
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FIG. 9: Magnetic field at the surface (see definition in the 
text) of a clean SC just before the entry of the first vortex 
for current approaching I c \. At d > A this magnetic field 
tends to the bulk critical field H c %. As d becomes smaller 
the SC can sustain larger magnetic fields in the vortex-free 
(Meissner) state. 



We turn to an estimation of the threshold current Is at 
which the barrier disappears. A more precise calculation 
would involve the Ginzbur-Landau theory. We follow the 
above analysis of H$£ Since the London theory is appli- 
cable at distances > ^ we estimate Is using 



dF , 



= 0. 



(13) 



We find numerically that at xq <C A the closed contour 7 
is well approximated by a circle with radius xo centered 
at x = z = 0. In the limit xq ~ £ -C A we can evaluate the 
functional F analytically as function of Xq. In Eq. (fTT|) 
for F v we can set exp(— \r — f\/X) — > 1, hence 



F v (xq) 



(/>IXQ 

327r/ioA 2 jq 



x 9(2x 



d6i I d02- — — — 7T-r- x 



sin ■ 



0- 



(14) 



Compared to F v , the stray term is negligibly small, 
F s ~ (^-) . The interaction energy with the ex- 
ternal current reads 



F ex t(x ) = - 
/to = 



2A 2 

X) 

dye' 



f(d/X), 



^(V^TT-y). (15) 



Using these formulas for and F ext we obtain from 
Eq. (fT3"|) the estimate 



(16) 



8^f(d/X) 



The dependence of /g on <i/A is hidden in the function 
fix), with /(a; -> 00) -► x"\ f(x 0) -> log(a:- 1 / 2 )+c, 



where 



0.3. For (1 > A we have Is 



<Pod 



In this 



case the magnetic field due to the external current at 
x = z = is [i?o(^ = 0)] 2 — > It is of the order of 
the second critical field H c2 . In the other limit rf<Awe 



have Is 



4>o 



4 Mo £log[A/d] 

£ c? <C A. In this regime we have Is 
Fig. (fTU)l we plot the phase diagram. 
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Note that this behavior holds for 
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FIG. 10: Phase diagram. For I < I c i vortices are thermody- 
namically unfavorable. For / > 7 c i curved flux lines become 
favorable, but an energy barrier opposes their entry until the 
current exceeds Is- A metastable phase exists in a narrow 
strip below the boundary J = I c i, shown in Fig. (0. 



B. Finite wire cross section 

Experimentally the wire carrying the external current 
has a finite cross section. Therefore it is important to 
include this effect in our calculations. Consider the rect- 
angular cross section as shown in Fig. (fTTj) . and assume 
current / flows uniformly in this cross section. We can 
write the external current as 



Jcxt 



X w Z<u 



d+- 



dd 



dz5(x — d)S(z — z)y. 



(17 J 

The modification to the magnetic field H = Ho + H v + H s 
occurs only in the first term, 



1 



^w^w J d 



dd 



dz 



H (x,y,z- z) 



i^d 



where Ho(x,y, z) is given in Eq. (j4|). Next we focus on 
the modification of the vortex dependent part of the free 
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FIG. 11: Rectangular wire cross section. 



energy F — F v + F s + F ext . Only the term F ext is modi- 
fied. Using Eq. (jTTJ) it is easy to find that 

F ext -> / dr z [°° dke- kd cos(kr z )(l-e- T ^) x 

1 ir io 



x 




Let us first specialize to the case of square cross section 
where the wire touches the SC, x w — z w = 2d, and com- 
pare this with a point like cross section x w = z w — ► 
which we considered until now (We ignore any electron 
or Cooper pair tunneling between the SC and wire). First 
we repeated the calculation of J c o and I c \. The results 
are roughly the same for both cross sections for d < A, 
and deviations up to 10% are obtained for d A up to 
d = 100A. In Fig. (fT2|) we compare the contours at I c \ 
as function of d for the two cross sections. We can see 
that zq changes by a factor of up to 1.6 for d < 90A. The 
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FIG. 13: Magnetic field at the surface (see definition in the 
text) just before the entry of the first vortex at / — *• J c i, for 
either zero wire cross section (x x = z w = 0, stars) or finite 
cross section (x x = z w — 2d, squares). 



z w — > oo can be treated analytically since the external 
field Ho is uniform at all x > (— d + ^f-). In this limit 
the maximal magnetic field at the surface before vortices 
penetrate is H sur f ace (I c i) — > H c \. For finite z w we cal- 
culated H sur f ace {Ici ) numerically with the result plotted 
in Fig. (fT4|) . Accordingly z w should not be too large in 
order to obtain the effect discussed here including the 
enhancement of the surface field in the vortex free state 
for a disordered surface. 

^surface 



H c i 
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FIG. 14: Magnetic field at the surface at current approaching 
Id as function of z m , for x w — 2d, d = A/2. 



FIG. 12: Dependance on d of the lengths xo and zo, charac- 
terizing the vortex contour, for zero versus finite wire cross 
section [x x — z w — (stars) and x x — z w — 2d (squares)]. 



magnetic field at the surface just below I c i is compared 
for the two cross sections in Fig. (fTHjl. At d ^> A it ap- 
proaches H c i in both cases, while for d <C A it is larger 
for point like cross section by about 10%. 

We estimate the typical values of the threshold current 
I c i. For the regime of interest d ~ A, we have I& of the 
order of — 1,6 w 55 ", lA . For A = lum this corresponds 

to current density of ~ 1mA — /jm~ 2 . 

Next we consider the dependence on z w for z w > 
x w = 2d which can be experimentally relevant. The limit 



IV. CONCLUSIONS 

In this work we studied solutions of London theory in 
a geometry where an external mesoscopic current flows 
parallel to a surface of a SC. Only above a threshold 
current I c q there exist solutions with curved flux lines 
entering and leaving the SC at the surface. At a larger 
threshold current, I c \, these solutions become energeti- 
cally favorable, however an energy barrier separates them 
from a vortex free solution. At a third threshold current, 
Is, this barrier disappears. To determine the current 
at which vortices actually penetrate the sample one has 
to account for the degree of disordered of the surface. 
For strong surface disorder the vortex can penetrate at 
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/ = I c i despite the presence of the barrier, due to large 
local magnetic fields produced at impurity sites allowing 
for nucleation of vortices. On the other hand for clean 
surface the entrance of vortices occurs at I — Is- 

By calculating those currents using a numerical solu- 
tion of the problem we conclude that for strong surface 
disorder the present geometry allows to achieve locally 
larger magnetic fields in the vortex-free state, as com- 
pared to the case of homogeneous magnetic field, pro- 
vided that the wire thickness is of 0(A). This can be 
potentially relevant for experiments in high-temperature 
superconductors which typically have extremely low val- 
ues of H c i . We argued that the effect of enhancement of 
the magnetic field in the vortex free (Meissner) state be- 



comes more pronounced in strongly anisotropic supercon- 
ductors, which is particularly relevant for layered high- 
temperature superconductors. 

We would like to thank Jordan Baglo, Walter Hardy 
and Cedric Lin for stimulating discussions. This work 
was supported by NSERC (ES & I A) and ClfAR (IA). 



APPENDIX A: MIXED TERMS IN FREE 
ENERGY 

We shall prove the vanishing of crossed term in the 
energy Eq. (TTJ) between H and H s [see Eq. flU], 



F {HoHa ) = Mo J d 3 r(H ■ H s + 9(x)X 2 (V x H ) ■ (V x H s )) = 0. (Al) 



From Eq. (0} we have H = H' + H sQ where H' = 0(-x)(Hi,d + H-i.-d) = V x Ai and 



r Iy , (x + d) 2 + z 2 



Correspondingly, we have F( HoHb -) = F{h' h b ) + F( HsoHa )- Consider the term F^ H ^ = ^ S x<0 d 3 rH' ■ H s . We will 

use the vector identity (V x A) ■ B = A - (V x B) + V • (A x B), with A = Ai, B = H s , and the fact that V x H s = 0. 
Then the volume integral can be transformed to an integral on the surface x = _ . However this integral vanishes 
because A\{Q~ ,y, z) = 0, hence F^ h > qHs ^ = 0. Now let us consider the term F^ HaoHs ) and define H s — V x A s . For 
the integral in the region x < we use the above vector identity with A = A s , B = H s q, and for x > we use the 
vector identity with A = H s q, B = V x H s . Taking into account that H s and H s q satisfy the homogeneous equations 
we obtain 

f d 3 r(V x A s ) -H s0 = f dS(l; x H- Q ) x , 
f d 3 r(H s0 ■ H s + A 2 (V x H s0 ) • (V x H s ))= -A 2 f dS(H+ x (V x #+)) (A3) 

Jx>0 J 



Here J dS — J^ ^ dy dz, and — H(x — O*, y, z). For x > we can use A s = — A 2 V x H s , which follows from 
London equation for H s . Next we use the fact that, by construction, H~ = + . This allows us to express 
Aj = — A 2 V x (H+ + H+), and combine the two terms of Eq. (|A3|) as 

F iHsoHa) = -MoA 2 J dS [(V x x (H- - H+) + (V x H v )+ x H M ] ^ . (A4) 
Now we use explicit forms of these factors: (H~ — H+ ) = — ^(dt+z*) ' i^7o)y = ^, 

( VxH °>y- 2^X 2 l dx JW e k 2 r(k 2 ) ' 



d 2 k 



Plugging these expressions in Eq. (|A4|) . one can readily obtain Fiu s0 H a ) = (without performing any integration), 
completing the proof for Ffjj jj a ^ = 0. 
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